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We investigate the discovery potential of cubic kilometer neutrino observatories such as IceCube 
to set stringent limits on the forbidden decays n — > v e v e and rj — > v T v T . The signatures for these 
decays are cascade events resulting from the charged-current reactions of v e , u T , v e and v T on nuclei 
in such detectors. Background cascade events are mainly due to v e 's from atmospheric jj,, K + , and 
K§ decays and to a lesser extent from atmospheric neutral current interactions with nuclei. A 
direct upper limit for the branching ratio r\ — > v e , T v e ,T of 6.1 x 10~ 4 at 90% CL can be achieved. 



PACS numbers: 13. 30. Eg, 14.40.Aq, 14.60.St 

I. INTRODUCTION 

The observation of the decay 77 — > vv or ir° — > vv would 
imply new and interesting physics. The r\ as well as the 
7T° have zero spin and odd intrinsic parity, i.e. J p = CP, 
and thus conservation of momentum and angular momen- 
tum require that the decay products v and v possess the 
same helicity. This decay provides an ideal laboratory 
to search for the pseudoscalar (P) weak interaction, be- 
cause only the P interaction allows the selection rule for 
the + — > 0~ transition for nearly massless neutrinos and 
antincutrinos. Other exotic effects such as the presence 
of a right-handed weak current through the exchange of 
a Z R would also allow such decays. The information de- 
rived from 7T° — > vv and 77 — > vv are complementary be- 
cause the former is sensitive only to the isovector neutral- 
current (NC) interactions while the latter is sensitive to 
the isoscalar NC interactions [l|. Furthermore, n° de- 
cays involve only u- and d-quarks while 77 decays addi- 
tionally involve the s-quark and perhaps other heavier 
quarks. If the Z° couples to a massive neutrino with the 
standard weak-interaction strength, the branching ratio 
(BR) for 7T° — > v T v T and 77 —¥ v T v T have maximum values 
of 5.0 x lCT 10 and 1.3 x KT 11 0], respectively at the v T 
mass upper limit of m„ T = 18.2 MeV/c 2 [If. It is note- 
worthy that BRs of « 2 X 10~ 18 and sa 2 x 1CT 15 are 
allowed within the Standard Model (SM) for 7r° —5- vvj 
and 77 — > vvj, respectively Q. 



II. EXISTING LIMITS 

To date no exclusive limits have been set on 77 —y vv in 
any experiment. The Particle Data Group (PDG)[|[ re- 
ports an inclusive upper limit of T(?7 — > invisible) /T(r] —y 
77) < 1.65 x KT 3 from the BES-II CollaborationQ, cor- 
responding to an upper limit on the BR for 77 — > invisible 
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of 6.0 x 10~ 4 . The BES-II results are inclusive results ob- 
tained by using 58 x 10 6 J/tp — > (fin decays. 

Possible 77 — > invisible decay products could be Light 
Dark Matter (LDM) particles or light ncutralinos. These 
LDM particles may have an adequate relic density to ac- 
count for the non-baryonic mass of the universe. Our 
estimated IceCube limits will be complementary to BES- 
II limits since the SM neutrinos would be a component of 
BES-II reported inclusive measurements. Limits on ir° —> 
v a v a (a = v e ,v^,v T ) are more common. An experimen- 
tal upper limit, r(7r° —> v e v e )/T{-n° -> all) < 1.7 x 10 -6 
at 90% confidence level (CL), was set by Dorcnbosch et 
al. The LSND Collaboration [f| has set an upper 

limit on the BR for tt -> v^ of 1.6 x 10~ 6 at 90% CL. 
In the tau neutrino channel, Hoffman has set a limit of 
r(7T° -> v r v T )/T(7r° -> all) < 2.1 x 10~ 6 at 90% CL 0. 
An inclusive search for tt° —> vv using K + — > tt + tt° has 
set an upper limit of 2.7 x 10~ 7 at 90% CL 0, (sec PDG 
for detail). 

III. CALCULATIONS FOR 77 -> v e , T v e , T 

Because of their enormous mass, cubic kilometer neu- 
trino detectors such as IceCube, offer a new opportunity 
to search for such exotic decays with competitive results 
compared to those obtained from accelerator-based ex- 
periments. IceCube, presently near completion at the 
South Pole, will contain 4800 digital optical modules 
(DOM) mounted on 80, 1-km strings. The active tar- 
get consists of approximately 4.2 x 10 37 le O atoms and 
8.4 x 10 37 H atoms. We have performed calculations 
using the CORSIKA Extensive Air Shower simulation 
code, version 6.72, to estimate the number of 77 mesons 
produced in the atmosphere The primary nucleon 
energy spectrum was calculated based on a sum of the 
power law approximations for the elemental primary en- 
ergy spectra, 

I{E) = Y,A i ^ Ai {EA i )-^. (1) 

The parameters and ja for i = 1,...28 primary 
nuclei with mass number Ai were obtained from corre- 
sponding approximations of balloon and satellite data 
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FIG. 1: Atmospheric v e (v e ), and v^ipy) energy spectra for 
two zenith angles compared with those of Gaisser and Honda 
[Till shown with solid and dashed lines. 



[loj . The resulting nucleon energy spectrum of expres- 
sion (1) can then be written as; 



1(E) = (0.110 ±0.006)£ 



-2.74±0.02 



(2) 



in units of (m 2 ■ s ■ sr ■ TeV)' 1 . The simulation program 
was tested by comparing the simulated neutrino and anti- 
neutrino energy spectra for two zenith angles (9 = and 
9 = 60°) with the corresponding spectra of Gaisser and 
Honda PH. 

Shown in Fig. Q] are the sum of fluxes of and v tll 
above 200 GeV at zero degree (open circles) and at 60° 
(open squares). Also shown are the sum of fluxes of v e 
and v e above 200 GeV at zero degree (solid circles) and 
at 60° (solid squares). These calculations, as shown in 
Fig HI agree well with those of Gaisser and Honda fllj ] . 
The v^Vfj) rate in FigJTJis an order of magnitude larger 
than v e (v e ) rate because pions and kaons decay mostly 
to v^(T>^) and not to v e (v e ). The largest sources of u e (v e ) 
are K% decay i.e. K± -> n°e^u e , BR = 5.08%, and K° 3 
decay, i.e. K\ — > n ± e T i' e , BR = 40.55% and to a lesser 
extent /i ± decay. Due to its long lifetime most /x ± reach 
the ground before decaying. Because the atmospheric 
i/ e (p e ) background is significantly less than the atmo- 
spheric v^(D^) background, we concentrate on searches 
for i] — > u e (v e ) and r\ — > v T (v T ). Interactions of high 
energy atmospheric neutrinos in IceCube can be classi- 
fied as events with a long muon track or as cascade events 
with very localized energy deposited in the detector. The 
muon events are due to charged-current (CC) interac- 
tion while the cascade events are mainly from CC v e and 
to a lesser extent from NC v e and v^. Neutrino absorp- 
tion in the Earth was also taken into account using the 



FIG. 2: The expected energy spectrum of v e (v & ) from decays 
of atmospheric K and \i is shown by open circles. The solid 
line shows the contribution of v e (v e ) from charm decays QiJ. 
The v e (i> e ) energy spectra from possible decay modes it — ► 
v e (v e ) (solid circles) and r\ — > v e (v e ) (solid squares) are shown 
assuming a 100% BR. 
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FIG. 3: v a integral energy spectra (see Fig 2) in the units of 
(cm 2 ■ year)^ 1 . 



neutrino mean free path X u = 1/ ' (Nap(&)(t(E v )), where 
Na is the Avogadro's number, p(9) is the average density 
of the Earth in g/cm 3 [ijj for a neutrino traversing the 
Earth at angle 9 and a is the i/-nucleon cross section at 
neutrino ener gy E v using a parton distribution functions 
from CTEQ6 [lj]. 
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FIG. 4: GEANT output for the spectrum of Fig. 3. The 
trigger requirement of at least 8 PMTs has been applied. 



In Fig. [5] the expected sum of the v e and v e energy 
spectrum from the standard K and fi decay modes is 
presented (open circles). Fig. [2] also shows the prompt 
u e (9 e )s from charm decay. This contribution becomes 
significant by 10 TeV and dominates at high energies [l6|. 
Also shown are neutrino energy spectra for the possible 
decay modes ir° — > v e v e (solid circles) and 77 — > v e v e 
(solid squares), assuming 100% BR. The energy spectra 
of neutrinos from a ir° and 77 SM-forbiddcn decays have 
nearly the same shape as the primary nucleon spectrum 
(7 ~ —2.7) whereas the energy spectrum of neutrinos 
from K and /_t decays is significantly steeper (7 ~ —3.58). 
This is because both ir° and 77 with very short lifetimes 
(t^o = 8.4 x 10~ 17 s and t v « 5.0 x 10 _19 s) do not inter- 
act and lose energy before decaying, while charged pions 
and kaons with much longer lifetimes interact substan- 
tially with the atmosphere before decaying. The cor- 
responding integral spectra are shown in Fig. [3j The 
flux of background neutrinos from the lower hemisphere 
(cos 8 < 0) is practically equal to the flux from the up- 
per hemisphere because neutrino absorption in the Earth 
is approximately compensated by the larger atmosphere 
depth in the case of upward-going neutrinos originating 
from the northern hemisphere. 

The flux of neutrinos induced by ir° and 77 decays 
slightly depend on atmospheric depth in the range of 700- 
1000 g/cm 2 . The CORSIKA generated neutrinos from 7r° 
and 77 decays with a CTEQ6 parton distribution model 
and the corresponding cross sections for CC and NC were 
calculated according to the formalism employed by Reno 
[l3| . The expected rate of detectable v e events for the 
IceCube detector was then calculated using a GCALOR 
simulation MC program [l4| ■ 
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FIG. 5: Measurable upper limits (lines) for the BR of tv , 77 — s- 
i/eVe decays versus neutrino energy for 5 years operating live 
time of the IceCube detector. The shaded areas show the 
statistical uncertainties. 



IV. GCALOR CALCULATIONS AND LIMITS 

The 3-momenta of these events for e and r leptons 
were written to a file which was then read by GEANT 
and electron and r transport with their accompanying 
hadron, a proton in this case, were simulated. The in- 
teraction vertices were distributed uniformly throughout 
the detector volume. The GEANT Ccrenkov code to- 
gether with the input IceCube geometry with average 
PMT quantum efficiency as well as an ice model with 
appropriate absorption and scattering [l5T | simulated the 
hit PMTs. The resulting trigger efficiencies for number 
of hit PMTs > 8 are 0.75 and 0.70 for 77 -> v e . T v e . T and 
the atmospheric background due to K and (i decays, re- 
spectively. Fig. 2] shows the energy distributions of the 
possible signal and the background events. We also esti- 
mated the atmospheric background due to cascades from 
ffj, NC interaction with nuclei as well as CC interac- 
tions where the muon is not detected due to edge effects. 
These types of events have a trigger efficiency of 0.30 and 
contribute only at a 15% level and have been included in 
the atmospheric neutrino background. The atmospheric 
v T contribution to the background are small below 10 
TeV (l6| . Furthermore, above a few TeV prompt v e (DJs 
from charm decay must be taken into consideration [l6j. 
The amount of MC data shown corresponds to one week 
of data taking with an 80-string IceCube detector con- 
figuration. Fig. [5] shows the expected upper limits on 
BRs for 7T° — > v e v e and ?/ — > v e v e decays that could be 
obtained from 5 years of measurements by the IceCube 
detector. The computations were performed at a 90% 
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TABLE I: Uncertainties associated with primary flux and neu- 
trino production cross sections. 



E v (GeV) 


10 J 


10 a 


10 4 


Primary Flux 


0.10 


0.11 


0.14 




0.15 


0.15 


0.15 


Overall 


0.18 


0.19 


0.21 



CL using the expression, 



where lnfc =< ln(E/E v ) >~ 2.75 ± 0.05. Table U shows 
these uncertainties for the primary flux and the neutrino 
production cross sections using the reported values of 
reference p7| . 

Uncertainties in table U are the most conservative esti- 
mates that contribute to the limits on rj — > v e>T v^ T . As 
shown in Fig. [5l the most stringent limits are obtained 
from neutrinos above « 9 TeV. 



BR < 



y iBkgr + {^-sys) 

/„(7r°, i] — > vv) 



(3) 



Where Iskgr is the number of background events from 
K, 7r, /i and prompt charm decays. A sys is the system- 
atic uncertainty in the number of background events, 
and I w (TT°,n — > vv) is the estimated number of n de- 
cay events assuming 100% BR. Note in the above ex- 
pression and the figure, the systematic uncertainties are 
mainly from two sources, the primary cosmic ray flux and 
cross sections. These uncertainties are energy dependent 
and a full analysis of them would include uncertainties in 
flux calculations and uncertainties associated with parti- 
cle production cross sections. We estimated these uncer- 
tainties based on the energy-dependent flux uncertainties 
reported b 
Derome [if 



17| and those reported by 



Agrawal, ct al. 
Uncertainties for the neutrino production 
cross sections reported by the two references is 15%. The 
relative primary nuclcon flux systematic error is due to 
uncertainties in Equation 2 and Iskgr simulations and 
can be approximated by 



Al 

T 



O.l 2 + (0.021n( 



lOOOGeV 



0' 



(4) 



V. SUMMARY 



In summary, we have investigated the IceCube discov- 
ery potential for setting stringent limits on the n — > v e v e 
and n — > v T v T . Our studies show that a direct upper 
limit of 6.1 x 10~ 4 at 90% CL at neutrino energies above 
9 TeV for both n decay to two e neutrinos or two r neu- 
trinos can be obtained. This limit is complementary to 
the limit set by the inclusive n — > nothing measurements 
of reference 4. 
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